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The conventional cosmological perturbation tlieory has been performed under the assumption that we know 
the whole spatial region of the universe with infinite volume. This is, however, not the case in the actual obser- 
vations because observable portion of the universe is limited. To give a theoretical prediction to the observable 
fluctuations, gauge-invariant observables should be composed of the information in our local observable uni- 
verse with finite volume. From this point of view, we reexamine the primordial non-Gaussianity in single field 
models, focusing on the bispectrum in the squeezed limit. A conventional prediction states that the bispectrum 
in this limit is related to the power spectrum through the so-called consistency relation. However, it turns out 
that, if we adopt a genuine gauge invariant variable which is naturally composed purely of the information in our 
, local universe, the leading term for the bispectrum in the squeezed limit predicted by the consistency relation 

vanishes. 
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I. INTRODUCTION 



The measurements of the cosmological fluctuations have now become a crucial tool to probe the history of the early universe. 
To yield the theoretical prediction of the fluctuations, which are to be compared with the observations, the gauge-invariant 
O . perturbation theory is used. It is not, however, widely recognized that the gauge invariance in the whole universe with infinite 
volume should be distinguished from that in the local universe with finite volume. It is widely known that the gauge-invariant 
^ ' cosmological perturbation can be realized by completely fixing the coordinates. In the case of the gauge invariance in the whole 
universe, where we request the invariance under normalizable gauge transformations, the conventional gauge fixing conditions 
are sufficient to uniquely determine the coordinates. By contrast, if we require that the coordinates within our local observable 
universe should be completely fixed only from the local information contained therein, the conventional gauge conditions are 
not sufficient. To fix the gauge locally, these conditions must be supplemented with appropriate boundary conditions [ 1 ] . These 
boundary conditions are fixed by requiring the regularity at the spatial infinity in the conventional approach. We should remind 
that, in the actual observations such as measurements of the Cosmic Microwave Background (CMB), we can observe only 
^—1 , our local universe and we cannot consult the regularity at infinity. Therefore, the observable fluctuations should be calculated 
in the request of the invariance under both the normalizable and non-normalizable gauge transformations. We distinguish 
^ gauge-invariant quantities that can be constructed from our local observable universe, referring them as genuine gauge-invariant 
. variables. The observable fluctuations should be such a genuinely gauge-invariant variable. 
I ' The importance of the genuine gauge invariance is also highlighted in the context of the infrared (IR) divergence problem. The 
I adiabatic vacuum, which yields the scale-invariant spectrum, is supposed to be a natural vacuum in the inflationary universe. 
. ^ ■ However, once the interaction turns on, it is not manifest whether the adiabatic vacuum is stable or unstable against the IR 
, contributions in loop corrections. It is indeed known that the loop corrections calculated in the conventional perturbation theory 
H ■ yield the logarithmic divergences ll2l-l2lll. Our previous works ll22i l23ll indicate that this divergence is an unphysical artifact 
. . .' in the case of single field models of inflation. We found that the initial quantum states should satisfy certain conditions if 
we require that the loop correction to the two-point function of a genuine gauge invariant vaiiable should be free from the IR 
divergence. Our interpretation is that these conditions are the requirements that the initial states are invaiiant under the residual 
gauge transformation. Namely, as long as the initial quantum states respect the gauge invariance in our local universe, the IR 
contribution to the loop correction should be suppressed. Hence, we referred to these conditions on the initial quantum states as 
the gauge invariance conditions. 

In this paper we reexamine the primordial non-Gaussianity in single field models, taking into account the importance of gen- 
uine gauge invariance of the selected variables and the gauge invariance conditions on the initial quantum state. A conventional 
prediction states that the bispectrum in the squeezed limit is related to the amplitude and the spectral index of the power spec- 
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trum, which is often called the consistency relation. However, it will turn out that the bispectrum of a genuine gauge invariant 
variable, which is the counterpart of the conventional gauge invariant curvature perturbation, is more suppressed in the squeezed 
limit compared with the prediction widely accepted as the consistency relation. 

This paper is organized as follows. In Sec. |II] we review a method to introduce genuinely gauge-invariant variables and the 
gauge invariance conditions. The basic perturbation equations will be also provided there. In Sec.|lII] for illustrative purpose, 
we first briefly summarize the usual argument for the consistency relation, which uses the information about the region outside 
our observable universe. After that, we re-investigate the bispectrum in the squeezed limit, using a genuinely gauge-invariant 
variable and imposing the gauge invariance conditions on the initial quantum state. Section|IV]is devoted to the conclusion. 

II. GAUGE-INVARIANT PERTURBATIONS IN LOCAL UNIVERSE 

Gauge invariance is guaranteed if we compute after complete gauge fixing. However, complete gauge fixing is not so straight- 
forward as mentioned in the introduction. Let's first impose conventional comoving gauge conditions 

<5</) = 0, (2.1) 

and 

SY^ = = d,6Y, , (2.2) 

with 

h,^=e'iP+C)[,s,^^^ , (2.3) 

where we denote the scale factor by e^. The time coordinate is then fixed completely by the first condition (12.1b . The second 
conditions ( |2.2| i. however, do not fix the spatial coordinates completely if we just consider the local universe and adopt gauge 
conditions only in this region. Namely, there are still residual gauge degrees of freedom as long as only the local metric is 
concerned. If we could use the information stored in the region outside our observable universe to fix the coordinates, these 
residual gauge degrees of freedom would disappear In spite of this, once we accept the contamination from the outside of our 
observable universe to define gauge invariant variables, such quantities cannot correspond to actual observables in a strict sense. 
Furthermore, such contamination from outside of our observable universe may cause artificial infrared divergences. We therefore 
need to fix the gauge completely using only the quantities residing in our observable region, but this is technically complicated. 

In our previous paper we pursued an alternative way to introduce genuine gauge-invariant variables, focusing on the scalar 
under three-dimensional diffeomorphism, e.g. three-curvature on the time-constant hypersurface. Such a quantity should be a 
genuine gauge invariant if the spatial position is specified in a coordinate independent manner In order to specify a point in a 
coordinate independent manner, we used the geodesic normal coordinates spun from a representative observer There remains 
the ambiguity in choosing the position of this representative observer and his/her frame. This dependence, however, should 
disappear in the evaluation of n-point correlation functions as long as the initial quantum state respects the translational and 
rotational symmetries. 

We computed the two-point correlation function of thus defined genuine gauge-invariant operator at one loop level, focusing 
on the effect of IR modes. Although the IR contribution to the two-point function is divergent in general if the spectrum of the 
linear curvature perturbation is scale invariant or red, we found that, if we choose the initial quantum states appropriately, this 
divergences disappear. We derived the conditions for the absence of IR divergences and confirmed that such quantum states exist 
at least up to the second order in the slow roll expansion. 

Now it is profitable to ask the question how to interpret the divergences for generic quantum states. Since Fourier modes 
of ( basically become constant in the IR limit, they can be almost absorbed by the overall scale transformation of the spatial 
coordinates, which is a part of the residual gauge degrees of freedom mentioned above. Although it is just a gauge, the fluc- 
tuation of this mode causes infinitely large variance of the two-point function in conventional computations. Hence, it was 
speculated that this divergence occurs because generic initial quantum states do not respect the expected invariance under these 
non-normalizable gauge transformations. Taking these in mind, we refer to the conditions necessary to avoid IR divergences in 
n-point functions as th e gauge invariance conditions. In this section, we briefly review these gauge invariance conditions derived 
in our previous works 11221 l23ll . 

A. Gauge-invariant operator 

The construction of genuine gauge-invariant quantities can be achieved simply by making use of scalar quantities under three- 
dimensional spatial diffeomorphism. In our previous works ll22l l23ll . we picked up the scalar curvature "i? as such a scalar 
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quantity. Although scalar quantities also vary under the gauge transformation because of the change of their arguments x', this 
gauge ambiguity does not appear in the n-point functions of these quantities if the arguments are specified in a gauge-invariant 
manner In order to specify the arguments of the n-point functions in a gauge invariant manner, we measure distances from an 
arbitrary reference point to the n vertices by means of the geodesic distance obtained by solving the spatial three-dimensional 
geodesic equation: 



dA2 



dx^ dx'' 



(2.4) 



where ^T^jk is the Christoffel symbol with respect to the three dimensional spatial metric on a constant time hypersurface and 



A is the affine parameter We consider the three-dimensional geodesies whose affine parameter ranges from A 
initial "velocity" given by 



dx'(X,A) 



dA 



-C(A=0) 



A=0 



-57(A=0)/2 



to 1 with the 



(2.5) 



We identify a point in the geodesic normal coordinates X* with the end point of the geodesic, .t' (X, A = 1). 

During inflation, the scales relevant to the current observations go far outside of the inflationary horizon, which leads to 
the natural assumption that the observable scale is much larger than the horizon scale l/H. In Sec. |III] we will evaluate the 
fluctuation with the modes 1/iobs ^ ^ e''H, where iobs denotes the observable scale in the comoving coordinates. In this 
section, however, it is sufficient to restrict the expression for to the case with 1 / Lobs ^ k, as we are interested in the IR 
divergence. In this case T'jfc is suppressed since it contains a spatial differentiation. Then, the geodesic normal coordinates X' 
are approximately related to the global coordinates x' as 



x\X) =: X' + SxHX) 



(2.6) 



We use this simple expression in this section, and we defer its generalization to the case with 1 / iobs <C fc to Sec. IIIIBI 
Using the geodesic normal coordinates X\ a genuine gauge invariant variable can be constructed as 



mix) := m{t, x^{x)) = • • • d^Rit, x^)U.=x. 



(2.7) 



where t denotes the cosmological time and we introduced the abbreviated notation X := {t, X}. 



B. Solving non-linear perturbation 

In this paper we consider a single scalar field with the canonical kinetic term, whose action takes the form 



Ai f 



(2.8) 



where Mp\ is the Planck mass and the scalar field (f> was rescaled as — > 0/A/pi to be dimensionless. We use the ADM 
formalism, following Ref. ll24ll . The metric in the ADM form is given by 

ds^ = -N^dt^ + h,j{dx' + N'dt){dx^ + N^dt) . (2.9) 

The above action is then rewritten as 



M. 



d'^x , (2.10) 



where "R is the three-dimensional scalar curvature and Eij and E are defined by 



Eij = [dth,, - D,Nj - DjN,) , E = h'^ E,, . (2. 1 1) 



Here Di denotes the three-dimensional covariant derivative and the spatial indices are raised and lowered by hij . To calculate 
the non-linear corrections under the slow-roll approximation, it is convenient to temporally work in the flat gauge: 

k,=e^f[e'^] (2.12) 
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with the transverse-traceless conditions 



5f,^0 = d,5f^, (2.13) 

because all the interaction vertices are explicitly suppressed by the slow-roll parameters in this gauge ifHllSl- Here in this section 
we associate a tilde with the metric perturbations in the flat gauge to discriminate from those in the comoving gauge. The action 
in this gauge is given by 



S 



N- 



m=0 



dt<rx, (2.14) 



,IR,. 



where y^") := d'"y/d(?i". We use the symbol "S" as in |23tl to denote an equality which is valid when we neglect the terms 
that do not participate in the IR divergences. When we write down the Heisenberg operator in terms of the interaction picture 
field, this is equivalent to keep only terms without spatial and/or time derivatives and the terms that include only one interaction 
picture field with differentiation. 

We also adopt the slow -roll approximation for the background evolution. To characterize the deviation from the exact de Sitter 
evolution, we use the horizon flow functions: 



£0 



H 



-rn-\-l 



1 d£„, 
£,„ dp 



for m > 0, 



(2.15) 



where H is the Hubble parameter and Hj is its value at the initial time. The horizon flow functions are related to the conventional 
slow -roll parameters as shown in Ref. ll25ll . Hereafter, assuming that the horizon flow functions Em with m > 1 are all small of 
0(e), we neglect the terms of 0(e'^). 

As studied in Appendix of Ref. ||24}l . the curvature perturbation in the comoving gauge C, is related to the fluctuation of the 
dimensionless scalar field (divided by Mpi) in the flat gauge Lp as 



IR 1 

C ~ Cn + 2^2C + CndpCn , 



(2.16) 



where we have introduced '■= — (/9/0)<p as in Ref. Il24tl . In the calculation of the n-point functions, we solve the evolution 
equation (Heisenberg equation) for the operator (f iteratively to express ^p in terms of the interaction picture field for ip. Variation 
of the total action with respect to ip yields 









+ pj 





yim+l) 



Variations with respect to the lapse function and the shift vector, respectively, yield the Hamiltonian constraint: 



1/(1 



-p"" + 2pd,N' + <j)ip + -p- 



IR 



(2.17) 



(2.18) 



and the momentum constraints: 



IR 



These constraint equations are solved to give 



4p 



~ IR • 1 



Substituting Eqs. ( |2.20| l and ( |2.21t into Eq. ( 12.171 ). the evolution equation of is recast into a rather compact expression. 



IR, 



= -2p 



,-2p 



4 p'^ 



where, changing the time coordinate from t to p, the differential operator L is given by 



a2 + (3-£i + e2)ap- 



(2.19) 

(2.20) 
(2.21) 

I, 

(2.22) 
(2.23) 
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Expanding (n and 6jij as 

Cn=1p + Cna H , = + 5%,2 H , 

the equation of motion ( 12.221 ) is reduced to 

= , (2.24) 



IR 



3 . e"2p 



'CCn.2 » -£ie2V'5pV' - 7e2e3V''^ - 2£i'0^^T^9^'0 71^5^1' did -jxp . (2.25) 

The first-order equation (I2.24t gives the mode equation for the interaction picture field and the second-order equation (12.251 ) is 
integrated to give 

Cn.AX) « (y + 6) + ei'/'ap^ + + S2)i^dpij + (5C„,2 + AaVX^p " ^'^xOV- + if^Tl'^.^x. V' • (2.26) 
Here, (n,2 includes the non-local term: 

'3 



(2.27) 



We kept the terms with second derivatives on the right hand side of Eq. ( 12.25b because they are also necessary to keep the terms 
in Cn,2 which have only one interaction picture field ip with differentiation. It should be emphasized that the homogeneous 

IR, 

solutions ^2^'^ and X2ip{dp — X^dxi)'4', in the sense that they satisfy £ • • • « 0{e'^), can be added to (n.2, where the time 
dependent functions ^2 and A2 should be O(e^) and their derivatives should be O(e^). 

C. Gauge invariance conditions of initial state 

Now, we address the gauge-invariance conditions to be imposed on the initial state, where initial conditions have two aspects. 
One is how we choose the positive frequency function for the interaction picture field in solving the mode equation ( 12.24b . 
The other is how we construct the non-linear Heisenberg field iteratively from the interaction picture field, i.e. how we add 
homogeneous solutions in Eq. (12.26b . The interaction picture field ip is quantized as 

(^e--V^.(p), (2.28) 

with 

V'fc(p) = Vk{p)ak + vl{p)a^^^. , (2.29) 
where the creation and annihilation operators satisfy the commutation relation 



In Ref. 12311, it was shown that the regularity of the one-loop corrections to the two-point function implies that the positive 
frequency function should satisfy 

[(1 + ^l2) dp-{l + X2)X'dx^ + /ii - 2£^V3] Vk{p)e'''-'' = -(1 + A2)i?fcffe(p)e^*"^ , (2.30) 
where Dk and £fc are defined as 

and we introduced the time dependent functions p^ (i = 1, 2, 3) written in terms of the horizon-flow functions as 



Dfe fc • Sfc + 3/2, j:^:=dl + {3-ei+S2)dp + ^k\ (2.31) 



1 3 

Ml £1 + -£2 + 2^2 , Ai2 £1 (1 + £1 + £2) + A2 , P3 := ^£2(2£i + £2) ■ (2.32) 

The condition ( |2.3Qb should be satisfied for all wavelengths and restricts the behaviours of not only the IR modes but also the 
UV modes. This condition can be shown to be consistent with the mode equation ( 12.241 ). At the leading order in the slow-roll 
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approximation, this condition, together with the usual condition on the UV behaviour, completely determines the initial state to 
the Bunch-Davis vacuum. However, if we extend our argument to the higher order in the slow roll expansion, Eq. (12.30) gives a 
non-trivial condition on the initial state. It might be surprising that such a gauge-invariance condition (= IR regularity condition) 
significantly constrains the allowable initial state. 

A gauge invariance condition on the coefficients of the homogeneous solutions is again derived from the regularity conditions, 
although it is not sufficient to completely fix both ^2 and A2. Imposing this gauge invariance condition, the potentially divergent 
terms in loop corrections add up to be total derivative terms with respect to k and hence they cease to contribute to the IR 
divergences. As discussed in Ref. [23], the time-dependent functions ^2 and A2 should also satisfy another relation to keep the 
commutation relation between and its conjugate momentum satisfied. Assuming that ^2 and A2 are appropriately determined, 
we do not discuss this relation any further because their explicit forms are not important for our current discussion. 



in. PRIMORDIAL NON-GAUSSIANITIES 



In this section, we evaluate the primordial non-Gaussianity in the geodesic normal coordinates for the vacuum that satisfies 
the gauge-invariance conditions described in the preceding section. Our main focus is on the so-called consistency relation that 
the bispectrum in the squeezed limit is given by the amplitude and the spectral index of the power spectrum. We will show 
that, if we consider the genuine gauge-invariant correlation functions as defined in the preceding section, the leading term in the 
squeezed hmit that arises in the usual computation of the bispectrum vanishes, and the well-known consistency relation does not 
hold. 



A. Standard consistency relation 

Now, we quickly review the derivation of the usual consistency relations for tree-level bispectra (CfciCfc2Cfc3) and {S'y^_^Ck2Ck3) ■ 
We are interested in the squeezed limit ki k2 — ks, where how we treat long wavelength modes becomes most essential. In 
this case the scale corresponding to fci crosses the Hubble horizon much earlier than the others, so that (^.-^ and Sjki become 
approximately constant by the time when the scales corresponding to k2 and fca cross the horizon. Therefore the mode with 
fci affects the evolution of the other modes mostly through deformation of the background geometry. The main effect of this 
deformation is taken into account by shifting the wavenumbers as 



fca ka= e '''=1 



ka = ka- [ Cfci + 7;^lki ) fca H , (3.1) 



2 

where a = 2, 3. This modification influences the fluctuations (k^ and Cfc,, leading to the replacement 

Cfc. (det/i)^y^Cfc„ Cfc„ - Cfci (ka-dk^ + Cfc„ - ^ka-5jki -dkXk^^ , (3.2) 

where the factor (dct h)'^^^'^ takes care of the change of the volume due to the modes with fci. Using the relation ( |3.2| i, we can 
easily evaluate the bispectrum for C, in the squeezed limit as 

(CfciCfcaCfea) — ^Pki (Slogfca + C^log fes + 3) (Cfc^ Cfc,,) 

~ -Pfc,<5(3'(fc2 + fca) [aiogfe, + 3] Pfc, ~ -,5(3) _ i)p^^p^^ ^ (33) 

where Pk is the unnormalized amplitude of the two-point function for C,k defined by 

(CfeiCfc.) =^^'Hfci + fc2)Pfei , (3.4) 

and ng is the spectral index. Thus, one can see that the three-point function in the squeezed limit is represented by and Pk- 
Equation ( 13.3b is called the consistency relation in single-field models of inflation. 
In a similar way, expanding the gravitational waves as 

hi],k = ^elj{k)5Yk. (3.5) 

with the polarization tensors defined by 

e^.%(fc) = = fc*e?(fc), (3.6) 
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the three-point function for SjCC is given by 

\a=l 

where P^*^' is the ampHtude of the two-point function for the gravitational waves. These three-point functions ( I3.3l l and ( I3.7l i 
were first calculated by Maldacena in Ref. ll24ll and these results were extended to more general single field models in Ref. Ii26ll . 



fe2-e''(-fci)-fc2 



(3.7) 



B. Gauge-invariant curvature perturbation 



In the following, we would like to consider the tree-level three -point function for genuine gauge-invariant variables. To 
compare the results with the standard computation described in the preceding subsection, it would be better to use C instead of 
i? as a variable, but ^ itself does not transform as a scalar on a <-constant hypersurface. Therefore the curvature perturbation 
evaluated in the geodesic normal coordinates; 



X{X) := C(p, x\X)) = C{X) + Sx'd,CU,=x' + ' 



(3.8) 



is not a genuine gauge invariant variable. However, this quantity %{X) is related to the genuine gauge invariant variable 
truncated at the required order in calculating the tree-level bi-spectrum: 



as 



(3.9) 



(3.10) 



Here we neglected the terms quadratic in gravitational wave perturbation. For the modes with k 1/ Lobs, using the Fourier 
decomposition of the variables, spatial derivatives can be replaced with the wave number, and hence Eq. ( 13.10) implies that % 
is equivalently gauge-invariant as ^R. Therefore, in the succeeding subsection we consider the three-point function for % in 
Fourier space. 

Now, we generalize Eq. (12.6b to be consistent with our situation 1/Lobs ^ <C aH. In spite of this effort, it turns out later 
that the deviation from the formula (12. 6t does not contribute to the leading order terms in the squeezed limit. Below, we mainly 
consider Fourier components of three-point function. One may wonder if the use of Fourier components is in conflict with 
our philosophy that genuine gauge invariant quantities must be defined only in terms of the information contained in our local 
observable universe. Here the relevant modes are restricted to sufficiently short wavelength mode in the sense that fc ^ 1/ Lobs- 
In this case, we can take Fourier components after multiplying a smooth window function. 

Focusing on the contribution from a single mode with the wavenumber fc, the geodesic equation is approximated by 



dA2 



[2(fc • X)X' - X^k'] Cke^^' 



X 



s=± 



1 



(k ■ X)X^^\^ - -k'X^X^5-,,k.k 



(3.11) 



where we neglected the non-linear terms, that do not contribute to the tree-level bi-spectrum. This equation is integrated form 
A = to 1 with the initial condition ( 12.51 ). The result of the integration ( 13.1 lb would be given in the form 



x\X) 



e-«(^) e 



^yix)/2Kx^ + Axn^k,Sjk]{X) + SxliijkKX) + dxfr[Sjk]{X), 



where we introduced 

Ax'i^k, S^k]{X) := [MX) - MO)] X' + [S^k.jiX) - S^k^jiO)] X^ =: Ax{ + . 



(3.12) 



(3.13) 



The third and fourth terms in (13.12) are corrections due to the curvature perturbation and the gravitational waves contained in 
the Christoffel symbol, respectively. As long as the contribution from a single mode ijjk is concerned, 5x]^ is given by 



SxliMiX) = ~ii:k[2{k- X)X' - X^k'] [ d\ [ dA 

^0 Jo 
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= -V/e [2(fc • X)X' - X^k'] / dk' dk 



:Ji^ik"-X 



{L)i ik"X 

^ 1 



where k" is to be understood as k"k/k. The tensor contribution Ja;^ can be similarly given by 



5xlr[5ik\{X) = Y,5^i^ / dfc' / dfc 



s=± 



{k-du,,)e^',di,,~-k^e^^,di,,dl 



ik" -X. 



(3.14) 



(3.15) 



C. Primordial non-Gaussianity measured in the geodesic normal coordinates 

In the present and succeeding subsections, we calculate the primordial non-Gaussianity measured in the geodesic normal 
coordinates. First, we neglect A.t, 5xl and 6xt in Eq. ( 13.12b . In this case, using Eq. ( 12.161 ). the gauge-invariant curvature 
perturbation is given by 



IR 1 1 

« Cn + ^ezC' + Udp - X'^dx^)(:n - ^Sj^'X,d.C„ 



(3.16) 



In the succeeding subsection we will extend our discussions to incorporate the effects of Ax, Sxl and 6xt, where we will find 
that this extension does not alter our main conclusion. Substituting the expression of (n given in Eq. ( |2.26l l into Eq. ( 13. 16b , we 
obtain 



TR 1 



(3.17) 



up to second order in perturbation, where the time-dependent coefficients /i; with i — 1, 2, 3 are defined in Eqs. ( 12.32b . Now, it 
is obvious that the contributions from the gravitational waves are cancelled in the genuine gauge-invariant curvature perturbation 
% and therefore the contraction of % with S'^ij trivially vanishes. The second-order perturbation in gravitational waves (57y,2 
also contain the terms with two ?/'s that can yield non-vanishing terms in the three-point function for %%6^. However, both of 
ips in Sjij^2 are associated with derivatives. We therefore arrive at the conclusion that the leading terms in (CC^j) obtained in 
Eq. (13.7b do not remain, once we use % in stead of (. 

Next we calculate the bispectrum of %s. Introducing the Fourier modes of % as 



d^X 



.-'"■""Xip, X) , 



(27r)3/2 

we consider (^Cfci^Cfca^Cfeg) the leading order in perturbation. Expanding as = + K2 + 



. we have 



(3.18) 



(3.19) 



Here the last term in Eq. (13.19b does not contribute to the leading term in the squeezed limit, fci ^ k2, k^ because the power 
spectrum for the mode with ki does not appear. 

Using Eq. (13.17b . the first term on the right-hand side of Eq. (13. 19b is recast into 

(V'fc,V'fe/Cfc3.2> = (i'k.i^k, {(1 + A*2) {^dp^)k^ - (1 + A2) {^/jx^dx^i^)^^+ yl^'jj) - (^fc.V'fc.£fc3V3 (V'')fc3> , 

where we defined 



d^X 

(27r)3/2 



-ikX 



n 



(27r)3/2 



P2 



(3.20) 



with O — I, dp, X'^dxi to concisely denote the terms originating from (n,2- The interaction picture fields ipk-^ and ipk^ should 
be contracted with either V'pj or tpp^ contained in {tpOtp)]^- Using the commutation relation for Uk and aj,, the first term of 
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Eq. ( |3.19t is rewritten as 

(V'/=lV'fc2^Cfc3,2> 





f d^X ' 


1 (2^)3/2 J 


' (27r)3 _ 



5(3)(fc2 {(1 + ^i2)dp - (1 + \2)X'dx^ + Ml - 2£;V3} i'^^P-^ 



^_,(fc,+fc3).x^(3)(fc^ + pX {(1 + m)dp - (1 + A2)X'9x. } v;e'P-'' 



(3.21) 



where we defined K := ki + k2 + k^. Here, we inserted the last terms on the first and second lines, which cancel with each 
other 

Following a similar calculation, the second term of Eq. (13.191 ) is rewritten as 



r d^p 


r d^x ' 


1 (2^)3/2 j 


(27r)3 _ 



3-'^'''+''^->''5^'Hks + pX {(1 + f^2)dp - (1 + X2)X'dx^ + Ml - 2£p Va} ^pe'^'^ 



i{k2+k3)-X 



j(3) (fci + p)v,^ { (1 + ;,2)ap - (1 + X2)x'dx^ } t;;e^p- 



(3.22) 



Replacing X with —X, we find that this expression is identical to the complex conjugate of Eq. ( 13.211) with the exchange 
between k2 and k^, which means 



k2.2l 



(3.23) 



Here we noted that at superhorizon scales becomes constant and commutes with the time integrals in where i = 2, 3. 
It is, therefore, sufficient to calculate the first term of Eq. (13.191 1. Using the gauge-invariance condition ( 12.301 ), the first line in 
Eq. ( |3.21t can be recast into the simple expression: llsill 



(First line in Eq. ( 13.211 )) 



= -(1 + A2)WfciWfc2 



7 (2^)3/2 7 (27r)3 Ifca+Pji^plt'pe j 

-(1 + X2)\vk, ?vu,j (^^-5(3) (fc2 +p){p-dp + 3/2) {y;5(^\k, +k,-p) 

d^p 



J(3) (fe2 +p){p-dp + 3/2) 5(3) (fc^ + k3-p) 

+ 5(3) (fci + fc3 - p) {(p + fcl - fc2 + fes) • + 3/2} 5(3) (fc2 + p) 

(27r)-3/2(fci - fc2 + fca) -91^5(3) (K) 



h J J0J^ dp p 5(3) (p + fei + fc3)5(3) {p-k2) 



(3.24) 



In the third equality, we replaced p with — (p — fci + fe2 ^ fca ) in the half of the expression. The second term in the last line takes 
the form of total derivative and hence it does not contribute. The first term proportional to ki ■ OkS'-^^ (K) in Eq. ( 13.241 ) vanishes 
in the limit ki 0. The other terms on the same line are not suppressed at this moment. Combining the contribution from the 
second term of Eq. ( 13.191 ), which satisfies Eq. (13.23b . these terms, however, provide a factor {vk2V*i^-^-i-k3\ ^ "|fei+fe2|^fc3)' which 
again vanishes in the squeezed limit, ki 0. 

To evaluate the second line in Eq. ( |3.21t , we again note that v^^ commutes with for ki <C efH. Thus, we have 



(Second Une in Eq. dllB) = -2\vkA^VkA27T)-''/H^^\K) {£pj^+^^|/i3< - ^fc/M3<} 
which vanishes in the limit ki — )- 0. 



(3.25) 
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As for the last line in Eq. ( 13.21b . we notice that p can be replaced with — fci owing to the factor S^^'' (fci + p). Therefore the 
term with the p-derivative acting on v* is suppressed in the limit fci — > 0. Hence, the terms in the third line can be simplified to 

= -(1 + X,)\vkA^\vkJ'{27r)-^/^k, ■ OkS^'^K) . (3.26) 

This expression vanishes in the limit fci — > 0. Thus we find that all the leading contributions in the tree-level bispectrum 
Eq. (|3.19t vanish in the squeezed Umit. 



D. Remaining contributions from the transformation to geodesic normal coordinates 



Here we take into account the effects of Ax, Sxl and 5xt- We calculate the three-point functions, whose three momenta are 
given by p^, P2, and pg without specifying which one is fci. We first discuss the scalar contribution Aa;^ to the three-point 
function with cubic (s, given by 



d^X 



P2I 



., |2 



(27r)3/2 



"Pi I 



12 « /■ d^X 



-ip^ X f^-ip^ X _ ^-ip^ X 



01 f\yp2 ?P2 ■ dp, s^^^ (pi + P2 + P3) - S^^\P2 + Ps) 



(3.27) 



When Pi is identified with fci, the factor [(5'''^-'(pi +P2 +P3) — ^'^■'(P2 +P3)] vanishes in the limit fci — > 0. When P2 is 
identified with fci, there is a manifest suppression factor proportional to P2. When pg is identified with fci, the power spectrum 
\vki P is absent and Eq. (13.27b does not yield the dominant contribution in the squeezed limit. Therefore, for all cases, the 
contribution from Axj^ vanishes in the limit fci — > 0. A similar argument follows for the non-vanishing contributions of Axt 
to the three-point functions. 

Next, we consider the contributions from 6xl and 6xt- Again we evaluate a rather general expression: 



d^X 
(2^ 



d^X 



^^e~VrX±iL)r ,k"-Xg -^p,■X 

(27r)3 -Pi''= 

|2| 12^2,1 
^'Pil M -T 
Pi 



r dfc' / dfc" [2(pi • du")di„~f^di„] s^'\p,+p, - fc"), 

Jo Jo 



(3.28) 



where fc" should be understood as — fc"pi/pi. Suppression is manifest in the case with P2 = fci or P3 = fci. In the former case 
there exists a manifest suppression factor proportional to P2, while in the later case there is no contribution from the spectrum 
\vki P as before. We therefore examine the case in which Pi is identified with fci. In this case, it will be convenient to rewrite 
the last expression as 



d^X 
(27r)3 



e-'P--''{^j;p^6x'){^p^d,4>) 



' Pi 



k' 



2{p,-dp^)dl,^~p\{dp^_-dpJ / dfc' / dk"S^^Hp,+p,-k"). (3.29) 



Denoting P2 + P3 by p, the part related to Pi is then simply given in the form 

|2 rPi r^' \n, |2 



^g(Pi) Tdfc' / dfc"5(3)(p_fc") = ^g(pi) rW / dk"\5{p-k")5^^\p + p{) 
Pi ' Jo Jo Pi Jo^ Jo P 

giPl)S^^\p + Pi) : 



2(pi -pY 6{pi ~p) 



PiP^ 



Pi-P 



(3.30) 



where p := p/p, Pi := Pi/pi, and g{p) :— 2{p- dp^)dp^ — p\dp^- dp,). If we integrate the above expression for a small 
regional < e = 0{p) multiplying G'(pi), which is an arbitrary smooth function of Pi, it becomes 0(e|uppG(p)). The obtained 
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expression is suppressed by a factor e compared with the case when we consider an integral for the simple delta function, 
/ d'^PiG{pi)\vp_^ (p + p^). Since this argument does not depend on the details of the function g{p), a similar suppression 
can be easily found also in the contribution from Sxt- Now, we can understand that all the neglected terms in the preceding 
subsection also vanish in the squeezed limit fci — > 0. 

E. Discussions 

The suppression in the squeezed hmit does not occur in the usual computation of the bispectrum for What we have found 
is that the leading term in the squeezed limit, exhibited in the consistency relation (13.3b . is significantly modified if we measure 
the gauge invariant curvature perturbation by using the geodesic normal coordinates. This can be understood from the fact that 
the non-Gaussianity described by the consistency relation is originating from the deformation of the local coordinates due to the 
IR modes that crossed the horizon much earlier than the other two modes. As shown in our previous works 1I22I l23ll and also in 
Refs. ll27|j29tl . this deformation can be absorbed by the gauge transformation in the local universe. This leads to the conclusion 
that there should not be the net physical effect due to these IR modes in this case. (See also Ref. (2^, where a similar analysis 
is presented, while the authors seem to have the slightly different opinion about what are the actual observables.) 

If we were thinking of the case with fci <C 1/ iobs, this cancellation can be understood rather easily as mentioned above. 
By contrast, our result might be a little surprising because the longest wavelength mode fci is, in the present setup, supposed 
to be still within the size of our observable universe. Nevertheless, the cancellation seems to continue to work. In the actual 
observation what we are interested in will be the temperature fluctuation of cosmic microwave background, which is not exactly 
what we computed here. In the case of CMB anisotropy, observation is performed not on a constant time hypersurface but along 
the light cone. Furthermore the observable region is not well outside the horizon even at the last scattering surface. Despite of 
these differences, our result suggests that for single field inflation models the leading terms of bi-spectrum in the squeezed hmit 
may vanish also in the CMB anisotropy as long as we analyze the data by using observationally natural coordinates. 

IV. CONCLUSION 

In this paper we reexamined the so-called consistency relation, focusing on what are genuine observables in our local universe. 
If we knew the configuration of the whole universe, usual curvature perturbation ( has gauge invariant meaning. However, the 
region we can observe is limited. In this case C itself is not genuine gauge invariant. 

In our previous paper we proposed a genuine gauge invariant variable that can be locally constructed, e.g. three curvature 
on (^-constant hypersurface measured in the geodesic normal coordinates. We found that the two-point function of this quantity 
is free from IR divergence in single field inflation models only if the initial state satisfies what we called the gauge invariance 
conditions. Applying thus-obtained gauge invariance conditions, we investigated bispectrum in the squeezed limit, where one of 
the momenta of the arguments of the three-point function is much less than the others. 

In the squeezed limit, the bispectrum is known to be expressed by the power spectrum and the spectral index. Although 
we are not claiming that the result of this standard computation is wrong, this conclusion is derived by evaluating the three- 
point function for the so-called curvature perturbation, which is not a genuine gauge invariant variable in the sense mentioned 
above. We therefore re-investigated the behavior of the three-point function for an alternative quantity, which is genuinely gauge 
invariant. As a result, we found that the leading term in the squeezed limit vanishes, i.e. the bispectrum is more suppressed in 
the squeezed limit than we expect from the naive calculation based on the so-called gauge-invariant curvature perturbation. 

We addressed the local non-Gaussianity in single-field models with the standard kinetic term, but our argument will be easily 
extended to models with a non-standard kinetic term. Since the large local non-Gaussianity is reported in multi-field models, 
the extension to multi-field models are intriguing, too. For this extension, we need to keep in mind the different nature of the 
IR modes between in single-field models and in multi-field models ||30ll . In the latter case, IR modes include physical degrees 
of freedom as well as unphysical gauge degrees of freedom. Therefore such perfect cancellation as we observed in single field 
models will not be expected in a general multi-field case. We will report these issues in our forthcoming publication. 
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